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Abstract 

We study, for the first time, the phase structure of the Gross-Neveu model with a combi- 
nation of a (constant) gravitational and a magnetic field. This has been made possible by our 
finding of an exact solution to the problem, namely the effective potential for the composite 
fermions. Then, from the corresponding implicit equation the phase diagram for the dynam- 
ical fermion mass is calculated numerically for some values of the magnetic field. For a small 
magnetic field the phase diagram hints to the possibility of a second order phase transition 
at some critical curvature. With growing magnetic field only the phase with broken chiral 
symmetry survives, because the magnetic field prevents the decay of the chiral condensate. 
This result is bound to have important consequences in early universe cosmology. 
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Because of its remarkable properties of being a renormalizable and asymptotically free 
theory, the Gross-Neveu (GN) model has been often considered as a reliable scenario 
for studying basic phenomena in particle physics, as chiral symmetry breaking and the 
formation of composite bound states. The study of the model under external conditions — 
as non-zero temperature, non-zero fermionic number density, or an electromagnetic field — 
has been carried out in Refs. (see also the references therein). In particular, an 

investigation of space-dependent configurations at non-zero temperature has been done in 
Refs. M. The phase structure of the GN model in an external gravitational field has been 



studied in Refs. |T0|, O], and a discussion of the chiral condensate in two-dimensional curved 



space can be found in Ref. |T^. For models of the early universe, the relevance of the 
study that we are going to present here originates in the fact that it is very probable that 
primordial magnetic fields have existed in combination with a strong curvature. Recently, 
there have been intense discussions on the role of such magnetic fields in relation with 



infiationary universe models |13]. Thence the interest in considering the combined effect of 
both gravitational and magnetic fields, which can be extended to quite different situations. 
In the case of the early universe this combined effect can produce a significant increase in 
the number of created particles | ]14|| . 

As a concrete example, we investigate in this letter the phase structure of the GN model 
on the two-dimensional Euclidean sphere S"^ with a magnetic monopole inside (in the three- 
dimensional language), that is, we consider the combined effect in the GN model of a mag- 
netic and a constant curvature gravitational field. As we will see, the magnetic monopole 
manages to prevent the decay of the chiral condensate. In fact, when the intensity of the 
magnetic field grows, the critical curvature tends rapidly to zero, and chiral symmetry is 
always broken. Our starting point is the action of the GN model with the auxiliary field 
a = Xifjip /N on the sphere S"^ of radius r (curvature i? = 2/r^) 

S = j (^r^^ij - \a^^ + y a^) , (1) 

where if) is an A^-component spinor. After performing the 1/A^-expansion and calculating 
as usually the effective potential for the non-zero condensate (a) = const 7^ 0, we can 
interprete the appearence of the condensate as a process of dynamical mass generation, so 
that M = A(cr). To first order in the 1/A^-expansion, we have 

V{M) = ^ - y Tr In (7^V^ - M) , (2) 
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where V is the volume and Tr includes a summation over the spinor indices and the trace 
of the covariant operator modes too. For the calculation of @ we will use zeta function 
regularization. It is convenient to work with the derivative of the potential 



^1 . It., (rv. - M)-' ^ ^1 - ^^IV fM^ + ^ v^y (3) 



Calculation of the trace yields 

where the well-known expressions for the spinor spectrum on 5"^ have been used (eigenvalues 
A; = {l + lf with de generacies di = 2 (/+!)). With the notation Mq = M|^^q and making use 



of the zeta-function regularization procedure (see Ref. [|T5[ for recent review of the method), 
we obtain 



dV 



= ^ + —Reij{l + irM), (5) 
A^ vr 

reg 



dM 

where ip{z) is here the digamma function. Here standard zeta-function regularization has 
been employed, in the sense of converting the divergent /-sum into a derivative of an Epstein 
zeta function of the kind F{z] a, b) = I].^o[('^ + '^)^ + ^]~^ with respect to a at a = 1/2, being 
b = r'^M'^ + 1/4. After that one just makes the analytic continuation of this zeta function 
to z = 1, taking due care of the poles (see W^)- Note that, written in this form, the result 



coincides with the one obtained in Refs. [^, where a full explanation of the regularization 
procedure, as applied to the present situation, is given. 

Having in mind the subsequent extension of the procedure to include a magnetic field 
(which is our goal here) we will introduce — already in the present case — the following 

0. Doing so, for the renormalized potential we 



renormalization condition: ^ 



M=Mo,r-*oo 



get 



dV 



— [Re^(l + zrM) -ln(rMo)]. (6) 

TT 



dM 

ren 

A graphic representation of the dependence of M on r, which is implicitly contained in this 
expression — equated to zero — is given in Fig. 1 (case k = 0, i.e., the curve below). As has 
been discussed in [|ll| in detail, this curve shows the behavior that corresponds to a second 
order phase transition. The critical curvature is defined by the condition In(rcMo) = Re 
ip{^), and has the value = 0.42. One should observe that we are not talking here about 
an absolutely real phase transition (since the only difference affects the discrete symmetry 
ip — > 75'?/', cr — o")- A delicate point of the approach is that of the possible influence of 
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space-dependent configurations, which could modify the whole picture substantially, leading 
to the single phase where M 7^ necessarily. 

On the other hand, one could also think in modifying the method through the intro- 
duction of some explicit symmetry breaking phenomenon (as an external magnetic field, 
for instance), which effectively supports the phase with M 7^ 0, even in the absence of a 
gravitational field. We shall now investigate this possibility in detail, taking as illustrative 
example the case of a monopole placed in the center of the spherical manifold S*^ (in the 
three-dimensional point of view). This issue of combining a gravitational and a magnetic 
field over a compact manifold had never been discussed before. We will see in the follow- 
ing that a precise model for such situation exists, and that sensible and interesting results 
can arise from this case. Since it is obviously impossible to introduce a magnetic field in a 
two-dimensional curved space of constant curvature, we came to the alternative idea that 
consists in describing the situation from a three-dimensional viewpoint, e.g., embedding S*^ 
in R^. From the point of view of this embedding only the normal component of the three- 
dimensional magnetic field B will influence the dynamics of the problem. However, one 
cannot impose that |B| = const, on the sphere (a compact manifold) since the total flow of 
B through the sphere is equal to zero (this was the original problem). Nevertheless, we can 
obtain an homogeneous flow of the magnetic field through the sphere by introducing in it a 
magnetic monopole (since then div B 7^ 0). If the monopole is located in the center of the 
sphere, its magnetic field has the form: B = k^, k = const > 0, and it can be considered as 
the natural analogue of a constant magnetic field B distributed on the sphere. Actually, 
this trick of introducing a monopole as a mean to study effects of a homogeneous magnetic 
field has been used extensively to simulate the quantum Hall effect (see, for instance, Ref. 
Ill7| ). In our situation, the role of B is played by the total flow = Airk (or by the quantity 
k itself). We will also assume that Dirac's quantization condition is fulfilled, and hence 
k = 1/2, 1, 3/2, 2, . . . For simplicity, we shall take e = 1. 

The change in Eq. (0) due to the presence of the magnetic monopole comes from an 
additional contribution in the connection [V^, V,y], which yields a term of the form e^yk/r"^ 
and, as a result, 7^7'^V^V,y = — -| — I75-R. Now the operator contains the potential 

and coincides with the Hamiltonian for a spinorial charged particle in the field of a Dirac 
monopole |T8|. The monopole does not break the spherical symmetry of the problem. The 
eigenvalues will actually change, but in a simple way. The effective potential (more exactly, 
its derivative) for the GN model on a sphere 5*^ with a magnetic monopole in its center (as 
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described above), is given by 



dV 



dM 



A2 V M2 - V2 + f + 

reg 4 2 



M 2(/ + l) 

where s = ±1 corresponds to the two different chirahties and k describes the magnetic 



field effect. Performing again a regularization via the zeta function [|15l, and using the 



renormahzation condition ^ 



= 0, we obtain the renormahzed potential. 

Af = Mo , r — > oo , fc /i xed 



The condition ^ = leads to an implicit dependence of the dynamically generated mass 
M in terms of r and k, M = M{r, k) (in place of @: 

J2 + l + isVr^M^ + ks' - F) = 41n(rMo). (8) 

s,s'=±l 

In Fig. 1 the dependence of y = M/Mq on x = tMq is represented for the first three values 
of fc, namely = 0, 1/2, 1 (lower to upper curve, respectively). The corresponding curve for 
A; = 5 is depicted in Fig. 2. Here, the loss of monotonicity through the appearence of a pick 
for higher values of k (starting at = 1) is clearly confirmed. Fig. 3 shows the value of 
the critical curvature Xc as a function of k. As we see, with growing k the critical curvature 
decreases quickly, what leaves the phase with M 7^ only. In other words, chiral symmetry 
is always broken then. We thus observe that the magnetic field has a stronger effect on the 
phase transition pattern than the gravitational field. For i? — one can easily find the 
asymptotic form of M = M{R, k) from (|]) analytically (for fixed k) 

12Mi 36QM^ ^ V ; K J 

From this last expression we see explicitly how the magnetic field acts against the decay of 
the chiral condensate. However, when k is not too big there exists a finite critical value Xc 
for which M{xc) = 0. The phase diagram shows a phase transition of second order. 

Up to now k has been kept fixed. However, in the situation when i? — 0, in order to 
have a magnetic field, B, we must scale k as k = Br^, B = const, r 00. In this limit B 
is to be identified with the usual homogeneous magnetic field on a plane. Taking this limit 
[B/Mq = 0{1)) one sees that M(r, k) tends to a finite value M{B), defined implicitly by 
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This expression confirms a well-known result for the dynamical fermion mass in the GN 
model with an homogeneous magnetic field. The first curvature corrections to the dynamical 
fermion mass are as follows (here M{B) > Mq, A; ~ Br^, R < B,M^, and B is finite): 



1 2B^ + 



(11) 



Note that when 5 ^ 0, we have from (^) that ~ Mq^ + 

Summing up, after being able to introduce a constant magnetic field on the surface of 
a two-dimensional sphere (serving as a generic example of a compact, curved manifold), 
we have studied the phase structure of the GN model under the combined influence of a 
gravitational and a magnetic field. The critical curvature has been shown to be a quickly 
decreasing function of k. With growing intensity of the magnetic field, chiral symmetry 
is always broken and hence there is no massless phase. We conclude, as a consequence, 
that the role of the magnetic field in this problem is more significant than the role of the 
gravitational field. It would be of interest to consider the generalization of the above problem 
to four dimensions, e.g., to study D = 4 four- fermion models influenced by a combination of 
a gravitational field, of Friedman-Robertson- Walker type, and a magnetic field, as the one 
created by a primordial monopole. Under the perspective of such more realistic situation, 
the neat (albeit preliminary) results that we have obtained here acquire a deep cosmological 
importance. 
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Figure captions 



Fig. 1. Curves showing the dependence of y = M/Mq on x = tMq represented for the first 
three values of k, namely k — 0, 1/2, 1 (lower to upper curve, respectively). 

Fig. 2. The corresponding curve, as in Fig. 1, for k — 5. Loss of monotonicity caused by 
the appearence of a pick for higher values of k (starting at A; = 1) is clearly observed. 

Fig. 3. The critical curvature Xc represented as a function of k for the values compatible 
with Dirac's quantization condition. With growing k the critical curvature decreases very 
quickly. This leaves the phase with M ^ only. In other words, chiral symmetry is always 
broken for k large. 
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Fig. 2: k=5 
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Fig. 3: x_c as a function of k 




